Dynamical central peak and spinon deconfinement in frustrated spin chains 
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Studying the dynamical spin structure factor in frustrated spin chains with spontaneously dimerized ground 
state we show that besides the gapped spin-wave excitations there appears at finite temperatures also a sharp 
central peak. The latter can be attributed to deconfined spinous, accounted well within the variational approach. 
The central peak remains well pronounced within the local spin dynamics and may be relevant for experiments 
on materials with ID frustrated spin chains. 
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Frustrated spin systems have been intensively investigated 
both theoretically an experimentally in last decades, offering 
novel phenomena and challenges as well as a broader view on 
strongly correlated electron systems. Among ID models the 
spin- 1/2 antiferromagnetic (AFM) Heisenberg chain (nearest- 
neighbor interaction J > 0) frustrated with the second neigh- 
bor AFM interaction J' > iH has attracted wide attention 
also due to its relevance to the quasi- ID material CuGeOa 
exhibiting the spin-Peierls transition atTsp — 14 K. For par- 
ticular parameters J' /J = a = 0.5 the exact ground state 
(g.s.) found by Majumdar and Ghosh (MG) f?) is doubly de- 
generate and dimerized with a spin gap to excited states. Such 
a spin-liquid state without a long range magnetic order has 
been shown to extend in a wider range around this point, i.e. 
in the regime a > ttc ~ 0.241. [3.-£5J Excited states at the 
MG point have been determined analytically and can be 
represented to a good approximation as a pairs of 5 = 1/2 
solitons or spinons with a gapped dispersion, the concept con- 
firmed by detailed numerical studies using the density-matrix 
renormalization-group (DMRG) method. [8] 

The dynamical properties of the frustrated J- J' spin chain 
have been so far mostly studied via the dynamical spin struc- 
ture factor S{q, uj) motivated again by the inelastic neutron 
scattering (INS) results on CuGeOa.lS] At T = the con- 
tinuum of 5 = 1 excitation in S{q, ui) in a wide range of a 
can be well represented with a pair of spinons, [8, 10] in par- 
ticular if the phenomenologically introduced matrix elements 
are taken into account 1 1 1 ] in analogy to the basic a ~ Q 
AFM Heisenberg model. 1 12] So far, there are very few the- 
oretical results on finite-T properties of frustrated system. It 
has been shown that the upper boundary of spinon continuum 
in S{q, uj) persists even at T > 0,(1^ while the maximum in 
the static structure factor S{q) exhibit a shift to incommensu- 
rate (7 < TT at larger a andT. 

In the following we present evidence that T > dy- 
namics of frustrated spin-chain model with the spontaneously 
dimerized g.s. exhibits several striking and rather unexpected 
features. Most evident, numerically calculated S{q,uj) re- 
veals at low but finite T > a sharp central w ~ peak 
well pronounced in the region g ~ tt, coexisting with the 
gapped two-spinon continuum known already from T = 
studies. [lofl The central peak has its manifestation in a 



unusual T-dependence of static susceptibility Xg^-n {T) with 
a maximum at T > 0. It shows up also in the local [q- 
integrated) S'l(w) as relevant, e.g., for the NMR spin-lattice 
relaxation. Using a variational presentation of excited two- 
spinon states and relevant matrix elements we show that the 
phenomenon can be directly traced back to spinons and their 
deconfined nature. 

In the following we study the frustrated spin model on a ID 
chain 



(1) 



where Sj are local S = 1/2 operators and the only rele- 
vant parameter is a = J' / J (we choose furtheron J = 1). 
The model has been invoked as the microscopic model for 
CuGeOs with a ^ 0.36 (realized above T < Tsp where 
lattice-deformation-induced dimerization is zero). But it as 
well represents the ID zig-zag spin system, example being 
double-chain compound SrCu02 within the opposite limit of 
large |a| ~ 5 - 10. d 

As the central quantity we calculate dynamical S{q,uj) at 
T > 0. We employ two numerical approaches. Finite-T 
Lanczos method (FTLM) ITtIi based on the Lanczos diago- 
nalization of small systems covers the whole T range but is 
restricted to system sizes < 28 whereby we use periodic 
boundary conditions (b.c). Finite-T dynamical extension of 
the DMRG (FTD-DMRG) method recently developed by the 
present authors 1 18] combines the DMRG optimization of ba- 
sis states with the FTLM method for dynamical correlations 
at r > and offers more powerful method for low T. The 
model, Eq.©, is here studied with open b.c. The reachable 
system sizes depend on T and the method shows good conver- 
gence, at least for low uj, for systems iV < 60 for T < 0.5 pre- 
sented here, with the typical subblock dimension m < 256. 
Concentrating on the low-w dynamical window the method 
is used as presented in Ref.[18], while high-w results are im- 
proved by the application of correction vectors increasing at 
the same time computation demand. The advantage of both 
methods is very good spectral resolution, so that typically only 
a minor additional uj dependent broadening of 5 ^ 0.02 at 
a; ^ and S 0.06 at higher uj is employed in presentations. 

Let us first present results for the MG model with a = 0.5. 
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While S{q, oj) at T = is rather well understood and inves- 
tigated numerically. lilQil we concentrate in Figs. 1 and 2 on 
T > FTD-DMRG results for different T/J < 0.5. The 
high-oj continuum appearing at T above the two-spinon 
gap Aq ~ 0.25 J 15J is qualitatively not changed from T = 
spectra. The evident new feature is the central peak at w ^ 
most pronounced atq ^ tt. Its width w is very narrow but still 
of intrinsic nature (being larger than the additional broadening 
6 ~ 0.02). It is evident that at fixed T the width w increases 
away from q = tt whereby the peak also looses the intensity. 
Still it remains well pronounced in wide region q > O.Ttt. 
The comparison of Figs. 1 and 2 also reveals that w increases 
as well with T and finally merges in a broader continuum for 
high T. 



be described well in terms of spinons as relevant excitations 
of the system and their deconfinement. At the MG point 
a = 0.5 the g.s. (for even N and periodic b.c.) has the energy 
Eq = —3NJ/8 and the wavefunction which can be written 
as the product of local singlets vj/g = [1, 2] [3, 4] • [iV - 1, A^]. 
It is doubly degenerate with the corresponding eigenstate 
having for one site shifted singlets. It has been akeady real- 
ized pi, l^l, ^ that lowest excitations can be well represented 
in terms of spinon states. In particular, the lowest branch of 
approximate triplet (S — 1,5^ = 1) eigenstates can be con- 
structed from the local triplet two-spinon states 

^b\p, m) -[1,2]... [2p - 3, 2p - 2] t2p-i [2p, 2p+l]... 

[2m - 2, 2m - 1] t2m [2m + 1, 2to + 2] . . . . (2) 




Figure 1 : (Color online) Dynamical spin structure factor S{q, u)) for 
the MG model at T/J = 0.1 within the whole range of g < tt, 
calculated by the FTD-DMRG method on a system of A*' = 60 sites. 
Correction vector improvement of high uj part was preformed only 
for q > 137r/15. 



where the first spinon is on site 2p — 1 and the second on site 
2m. Since the total momentum Q is conserved due to periodic 
b.c, the relevant two-spinon functions are 
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M 



^^(fc) = J_'^gi{Q+k)p+i(Q-k)mj,t 

p,m 



V'*(p,m), (3) 



where sums run over M — N/2 double cells. In the further 
analysis difficulties arise since ij/ {p, m) are not orthogonal 
even for distant |p — r7i| ^ 1 and furthermore forp ^ m. To 
find proper eigenfunctions we follow the procedure and nota- 
tion of Ref.[6J which for each momentum subspace Q yields 
nontrivial matrix elements 



9J2 1 

ITa Kk' + ^XQ{k,k'),(^) 



M 



where H = H - Eq, uj± = u}{{Q ± k)/2), uj{p) = (5/4 ■ 
cos 2p)J/2 are (approximate) single-spinon energies and 




(- + cosQcosfc)J. 



(5) 



Figure 2: (Color online) S{q — n,ijj) for different T/J = 
0.0, 0.1, 0.3, 0.5 where the broadening of central peak with increas- 
ing T is well pronounced. 

In the following we present the analysis showing that the 
emergence of the central peak in S{q,uj) at low T < J can 



The off-diagonal terms XQ{k,k'),hQ(k,k') (not presented 
here) emerging from nonorthogonality of ipgik) are the same 
as given in Ref.|6]. Within the triplet two-spinon basis, 
Eqs. (|2]i and (O, proper eigenstates (nevertheless not yet exact 
eigenstates of Eq. ([T])) are obtained via the diagonalization of 
Eqs.© and can be denoted 5'q(A:) (whereby k remains only a 
label and not a well defined wavevector). In contrast to i/jq (fc), 
\I'Q(fc) are ortho-normalized. The corresponding two-spinon 
excitation energies eQ(fc) are well approximated as the sum 
of two (deconfined) free spinons, i.e. egik) ^ ^gik), Eq.©. 

Our goal is, however, to understand low-T properties of 
S{q, u). Before discussing T > results, we first have to re- 
consider the T — spectrum So{q, w) which has been already 
interpreted in terms of two-spinon excitations. HE III Still 
the corresponding matrix element has been postulated so far 
only phenomenologically \ll. 23] in analogy with previous 
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works on the unfrustrated Heisenberg model. fl^] We note that 
at T = within the chosen subspace, Eq.(|2]i, we can express 



Soiq,u;) = IJ2 \{%(.k)\S+\^o)\'Siu; e,(fc)), 



and 



-iq 



2M 



E 

k 



(6) 



(7) 



Since Eq.(|7]i is an exact representation of the operator, 
we can evaluate matrix elements (q{k) — (5'* (fc)|S'+|4'o) by 
diagonalizing numerically equations, Eqs.(|4|i. Taking into ac- 
count that e*(fc) ^ £q{k) we can then evaluate 6*0(5, w) in 
the two-spinon approximation. Results within such a frame- 
work are presented for q = tt in Fig. 3 along with the full 
numerical results obtained via the T = FTD-DMRG (for 
T = identical to the more standard dynamical DMRG) 
evaluated within a system of = 100 sites. The agree- 
ment is very satisfactory except at the higher-w end where the 
obtained intensity is too low as well as Eq.(|4|i seem to gen- 
erate a high two-spinon anti-bound state (peak in So{tt,uj)) 
besides the free two-spinon dispersion, Eq.(|5]l. It should be 
noted that obtained Cq{k) is quite far from the oversimplified 
spinon picture with Cq{k) ^ l-[20] Still it is hard to find for 
it an appropriate analytical expression. JSi ll^l One possibil- 
ity is to neglect non-orthogonalities in Eqs.(|4]l which yields 
Cq{k) oc Corresponding "free" spinons results 

for So{q — TT, w) also presented in Fig. 3 show qualitatively 
reasonable trend (fall-off for higher w). Still they give an in- 
correct behavior at lower and higher cut-off due to divergent 
two-spinon density of states. 
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Figure 3: (Color online) T — dynamical spin structure factor 
So{q = 7r,t^) as calculated via the DMRG for iV = 100 sites (full 
line), numerically within the two-spinon approximation (dashed line) 
and using simplified ^q{k) (dotted line). 

The above agreement of numerical T = results with the 
description in terms of the two-spinon basis, Eq.(|2]l, gives firm 
support also to the interpretation of T > dynamics. In the 
low-T regime we are in S{q, lu) predominantly dealing with 



excitations increasing the number of spinons, — rig + 2, 
analogous to those in SQ{q, uj), Eqs. ^ and (|7]i. Their contri- 
bution analogous to T = Fig. 3 is evident also at T > in 
Figs. 1 and 2. 

However, in addition there are possible transitions between 
excited states conserving Ug. In particular, the matrix element 
lqQik,k!) = (*^+Q(fc)|5+|«'^(fc')) as introduced akeady 
in Ref.[201 is finite and nontrivial. Here ^Q(fc') are singlet 
two-spinons eigenstates. We evaluate 7gQ(fc, k') numerically 
assuming two-spinon approximation, Eq.©. Results show 
that elements are nearly diagonal, i.e., jqQ{k, k') ~ Sk'.k+Q 
leading in S{q, uj) to the contribution atco^ ^q+qi^ + q) ~ 
eg(fc). Since at low-T favored are lowest excited states, i.e., 
from Eq.© Q ^ 0,k ^ tt and Q ^ TT,k ^ with a 
Boltzmann weight p oc exp(— Aq/T) (where Aq ^ 6,^(0) = 
J/4). Numerical solution of two-spinon problem shows that 
^Q^k) ^ ^gik) ^ ^Q{k) consistent with the picture of un- 
bound (deconfined) spinons. Hence, the strongest transitions 
areata; ^ eq+Q{k + q) — eQ{k). This evidently leads atq ^ tt 
to a sharp central peak at a; with the strength increasing 
as oc exp(— Aq/T). 

From above perspective we therefore conclude that the pro- 
nounced central peak in Figs. 1,2 at q ^ tt confirm the pre- 
sented analysis of nearly-free or deconfined spinons as excited 
states at the MG point. On the other hand, even without the 
extensive calculations it is evident that the central peak can 
only appear if triplet and singlet spinon states are nearly de- 
generate again only possible for deconfined spinons. 

The emergence of the central peak is, however, not re- 
stricted to the MG point a = 0.5 but appears to be related 
closely to the existence of the spontaneous dimerization at 
a > ac and the spin gap Aq > 0. We tested numerically 
also the case a = 0.7 (only partly presented here) where the 
spm gap IS larger Aq - 0.4 J. (3\ Consequently also the cen- 
tral peak feature is even more pronounced and extended in the 
q space as well as persists to higher T. 

It is evident that the central peak has a substantial effect on 
the static susceptibility Xq{T), 



Xq{T) 



duj , 



(8) 



being sensitive to low-w dynamics. In Fig. 3 we show the 
FTLM results obtained on systems with = 28 sites for 
Xq{T) with various q and again a — 0.5. Most pronounced 
is the variation at g = vr where the g.s. value Xir(0) is deter- 
mined with the dimerization gap, i.e. X7r(0) cx I/Aq. Instead 
of naively expected monotonously decreasing Xq{T), we ob- 
serve in Fig. |4] simultaneously with the emergence of the cen- 
tral peak an increasing XTr{T) in the regime < T < T* 
whereby T* - Ao/2. On the other hand, for T > T* the 
fall-off is uniform with XTr{T) oc 1/T which is close to the 
characteristic critical q = tt behavior for the simple AFM 
Heisenberg model. [18?] Results for g < tt are quite analogous 
taking into account that the relevant spin gap is Aq > Aq. 

The effect of the central peak is visible also in local spin 
correlations — {i/N)J2nS{q,i^) as presented in 
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Figure 4: (Color online) Static susceptibility Xq [T) vs. T for differ- 
ent q < n for a = 0.5 as obtained with the FTLM. 

Fig. |5] They are accessible directly via FTD-DMFT by cal- 
culating local spin correlations {S^,S^)ui locating the site 
i ^ N/2 to avoid effects of open b.c. As well we can eval- 
uate them within the FTLM and periodic b.c. summing all 
S{q ^ 0, w) (g = contribution is delta function due to 
the conserved 8^^^). Clearly, in both approaches the diffusion 
contribution g is not represented correctly but it is ex- 
pected to be subdominant. i2lll FTD-DMRG results in Fig. 5 
presented for a = 0.5, 0.7 reveal a central peak at a; 
well separated from the higher-w two-spinons continuum as 
far as T < Aq. The peak gains the weight at T ~ T* and for 
T > T* steadily becomes broader, finally merging with the 
continuum for T > Aq. Moreover we observe for both a that 
S'L(aj = 0) is nearly constant in a broad range T* < T < 2 J. 
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Figure 5: (Color online) Local spin correlations Sl{oj) for a = 
0.5, 0.7 at T = 0.1, 0.3 as obtained with the FTD-DMRG method. 

It should be noted that within the simplest approximation 
(with q-independent form factor) the NMR or NQR spin- 
lattice relaxation should be a closely related to Sl{i^), i-e. 
the relaxation rate is given by 1/Ti oc ^^(i^ = 0). Fol- 
lowing above results we would obtain for considered systems 
l/Ti ^ const in a broad range T > T* similar to theoretical 
predictions for the ID (unfrustrated) AFM Heisenberg model 
and CuGeOa While the agreement for higher T > Aq 
with the Heisenberg model is not surprising the novel contri- 



bution of the central peak is that the validity of this universal- 
ity is extended to lower T > T*. It should be also reminded 
that such relaxation is far from the usual Korringa relaxation 
with l/(rTi) ~ const. 

In conclusion, we have shown that the frustrated spin chain 
as manifested within the ID J- J' model with a > ac reveals 
besides the gap in spin excitations at T = also very un- 
usual spin dynamics at finite but low T < Aq. The central 
peak which appears in S{q,uj) at q ^ n as well in the q- 
integrated local Sl{(^) is very sharp and dominates the low-o; 
response at low T. It is a direct consequence and the signa- 
ture of deconfinement of spinon excitations in such systems. 
It remains to be investigated whether such a behavior is re- 
stricted to the particular case of investigated model or there 
are other gapped spin systems with similar phenomena. As 
far as experimental relevance is concerned extensively inves- 
tigated CuGeOs above the spin-Peierls transition T > Tsp is 
interpreted with a frustrated spin-chain model with a ^ 0.36 
and could partly exhibit mentioned phenomena in spite of pre- 
sumably very small scale Aq < 0.02 J f?]. 
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